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Abstract 

The general analytic solution to the functional equation 



(x + y) 



is characterised. Up to the action of the symmetry group, this is described 
in terms of Weierstrass elliptic functions. We illustrate our theory by ap- 
plying it to the classical addition theorems of the Jacobi elliptic functions 
and the functional equations 



<h( x ) 


02 (y) 


03 0) 


03 (y) 


<t>4(x) 


04 (y) 


050) 


05 (y) 



0i (x + y) = 04(x)0 5 (y) + (/>4(y)05(aO 



and 

+ y) = ^2{x + y)<h(x)<h(v) + *3(z + y)Mx)Mv)- 



1 Introduction 



The purpose of this article is to describe the general analytic solution to the 
functional equation 

M x ) My) 
M x ) My) 



M x + y) 



(i) 



M x ) My) 
M x ) My) 

Although this equation appears to depend on five a priori unknown functions we 
shall show that ([!]) is invariant under a large group of symmetries Q and that 
each orbit has a solution of a particularly nice form, expressible in terms elliptic 
functions: 

Theorem 1 The general analytic solution to the functional equation ([][) is, up 
to a Q action given by ([7Z]-[7^) ; of the form 



M x ) 



Here 



M x ) 
M x ) 



and 



M x ) 
M x ) 



$0; i/ 2 ) 

$0; u 2 y 



x.1 \ a ( u ~ x ) CMx 



(2) 



a(z/)o"(a; 

where o~(x) = a(x\oj, a/) and ((x) = are the Weierstrass sigma and zeta 
functions. 

The group Q of symmetries of (|1|) will be described further below. Our proof is 
constructive and indeed yields more: 



Theorem 2 Let xq be a generic point for Then (for k = 1,2) we have 

= C(^k) - C(x) - C(u k - x) - \ k 



d y In 



M( x + x o) M(y + x o) 
M+i i x + x o) M+i {y + x o) 



y=0 



— A fe + y>, 

x / — ' 

1=0 



X 



l+l 



(/ + !)! 



and the Laurent expansion determines the parameters g\, g 2 (which are the same 
for both k = 1,2) characterising the elliptic functions of ([Jj by 



g 2 = - (F 2 + 6F 2 ) 



g 3 = 6F*-F? + ~F F 2 , 



and the parameters v k via F = —p(u k ). Further, we have 
M x + 2*o) = M2xo) e (A2 " Al)1 ' 

®{x; v 2 ) 



and 



<i>2k( X o) 4>2k(x ) 



<p2k{x + X ) 

hk+l(x + X ) J f( X ) \4>2 k +l( X o) 4>2k+l{Xo) ) \\ k 



v k ) 



Here the junction 



-\kX 



<t>2k( X o) 
^2fc+l( X o) 



4>2k(X()) 

hk+i(x ) 



<p2k{.X + Xo) <p2k(Xo) 
4>2k+l(x + Xq) 4>2k+l(x ) 



is in fact the same for k = 1,2. 



The term 'generic' will be defined below and we will give more expressions for 
the quantities appearing in the theorem. 

One merit of writing (p]) in this general form is that several different functional 
equations may now be seen as different points on a £?-orbit of ([I]). Thus for 
example 



M x + y) 
M x + y) 

A{x + y)[B(x)-B(y)} 



M x )Mv) 

(j>4(x)(p 5 (y) + fa(y)(j) 5 (x) 
A(x)A'(y)-A(y)A'(x). 



(3) 
(4) 
(5) 



are particular^ examples of ([I]). The functional equation for the exponential (|3|) 
corresponds to V\ = z/ 2 in our solution and the exponential comes wholly from 
Q. Particular cases of (H) have been studied in M and we shall determine (see 
Lemma 4) the general solution to (|j) as an application of our work. 

More interesting is equation (|5|) which has been studied by various authors 
with assumptions of even/oddness on the functions appearing |12|, [R| [L9| or 
assumptions on the nature of B [IB]. The general solution [|], ^ A(x) = v) 
now corresponds to to the limit z/ 2 — > together with a Q action. This will be 
illustrated later. 

Finally, when 4>\{x) = a(x), 4>2(x) = a(x)r(x), (p^x) = t(x), (p^x) = 4> 2 (x) 
and 4> 5 (x) = 4>' 4 (x) we obtain the functional equation 



a(x + y) = a(x)a(y) + T(x)T(y)ip(x + y). 

(j>l(x)(j) 5 (x), 



(6) 



1 These correspond to 

(a) <p 2 (x) = (f>i(x)4>4,(x) and <ft 3 (x) 

(b) fa (x) = 4>l{x) and (f> 3 (x) = 4>l{x), 

(c) 0i (x) = (j) 2 {x) = A(x), (f> 3 (x) = A'(x), <pi(x) = B(x) and <f> 5 (x) = 1. 
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The function ip(x) will be described in more detail in the sequel. This equation 
was studied by Bruschi and Calogero |3| and will be used in our analysis. 

Let us remark that both (H) and (0) may be viewed as limiting cases of the 
functional equation 



+ y) = V 2 (x + y)fa(x)fa(y) + V 3 (x + y)M x )Mv), ( 7 ) 

which a priori depends on seven unknown functions. Later we shall show how 
(0) may be used to solve this. 

It remains to place ([I]) in some form of context. The last decade has seen 
a remarkable confluence of ideas from completely integrable systems, geometry, 
field theory and functional equations that is still being assimilated. To make 
some of these matters concrete let us consider how such functional equations 
arise in the context of integrable systems of particles on the line. A pair of 
matrices L, M such that L = [L, M] is known as a Lax Pair; this is a zero 
curvature condition. Starting with an ansatz for the matrices L and M one seeks 
restrictions necessary to obtain equations of motion of some desired form. These 
restrictions typically involve the study of functional equations. The paradigm 



for this approach is the Calogero- Moser [11 system. Beginning with the ansatz 
(for n x n matrices) 

L jk = Pj5 jk + 9 {l-8 jk )A{q j -q k ), M jk = g [S jk ^ B {Qj—Qi)-0--Sjh)C{Qj—Qk)] 

one finds L = [L, M] yields the equations of motion for the Hamiltonian system 
(n > 3) 

H = — p 2 - + g 2 U(qj — q k ) U (x) = A(x)A(—x) + constant 

provided C(x) = —A'(x), and that A(x) and B(x) satisfy the functional equation 
(|5|). With this ansatz and assuming B(x) evenQ Calogero [Pj| found A(x) to be 
given by (0). In this case the corresponding potential is the Weierstrass jo- 
function: A(x)A(— x) = p{v) — p(x). The functional equation (|6|) is associated 
with a different ansatz and yields the relativistic Calogero-Moser systems H [21 



Similarly ([]]) arises from a more general ansatz M associated with equations of 
motion of the form 

q) = J^(a + bqj)(a + bq k )V jk {qj - q k ), 



This assumption can in fact be removed 
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which combines both relativistic (b ^ 0) and nonrelativistic (b = 0) systems 
together with potentials that can vary between particle pairs. This unifies, for 
example, Calogero-Moser and Toda systems |2T|, |22J. The relativistic exam- 
ples yield the functional equation (P while the nonrelativistic situation involves 
the functional equation 



M x + y) = M x + y)(M x ) - My)) + 



<j>' 2 {x) 



My) 
My) 



(8) 



The general analytic solution to (H) has yet to be determined although particular 
solutions are known. We remark that (HI), and after suitably symmetrizing (|8]), 
are particular cases of the functional equation 



N 



^M( x + y) 

i=0 



MM X ) M+i(y) 
MM x ) MMy) 



(9) 



with AT = 1 in the former case and iV = 2 in the latter. When 03i+2 = M+i 
Buchstaber and Krichever have discussed (|9|) in connection with functional equa- 
tions satisfied by Baker- Akhiezer functions [fTPf . 

Lax pairs are but one way in which functional equations are associated with 
integrable systems and we mention ]7|, 13), |9], [14| for others. There also appears 



a close connection between these functional equations and the elliptic genera 
associated with the string inspired Witten index |y| [17| . Krichever for example 
used the functional equation (Q) in his proof of the 'rigidity' property of elliptic 
genera [17] and it also appears when discussing rational and pole solutions of 



the KP and KdV equations |K| M. We feel this connection between functional 
equations and completely integrable systems is part of a broader and less well 
understood aspect of the subject that deserves further attention. 

An outline of the paper is as follows. First we will discuss the group of 
symmetries of (|l|). These will be used in the proof of theorem 1. Before turning 
to the proof we show in section 3 how the indicated solution does indeed satisfy 
(jl]), using this as a vehicle to recall some of the properties of elliptic functions 
that we will need throughout. Section 4 is devoted to the proof of theorem 1 and 
section 5 to that of theorem 2. Several applications of our theorems including 
the general analytic solution to and a discussion of (0) are then given in 
section 6. An appendix is given containing various elliptic function formulae we 
shall make use of. 

Various versions of Theorem 1 have appeared in unpublished preprints. In 
the form of 0i(x) only was stated. In we introduced the Q action to 
give Theorem 1 in its present form. In improving the proof of this we obtained 
Theorem 2, given here alongside the better proof of Theorem 1. 
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2 The Group of Symmetries 

We next describe the group Q of invariances of (|]). Theorem 1 gives a repre- 
sentative of each Q orbit on the solutions of ([!]) with a particularly nice form. 
First observe that a large group of symmetries Q act on the solutions of (Q). The 
transformation 

(10) 

clearly preserves ([!]) provided 



A + A' + A" = 0, U,V E GL 2 , and det C/ = c det V. (11) 
Further, (]!]) is also preserved by 




We will use these symmetries in our proof of theorem 1 to find a solution of (|1|) 
on each Q orbit with a particularly nice form. 



3 Illustration of the Solution 



Before proceeding to the proof it is instructive to see how the stated solution 
satisfies (|]). This will also allow us to introduce some elliptic function identities 
needed throughout. From the definition of the zeta function we have 



;in$(ac; u))'= -((v - x) - C(x) + C(v) 



(14) 



Thus 



= $(x; v)<fr(y; v) (ln$(y; z/))'-( In $(x; v))' 

= $(x; v)$(yi v) ((v - x) + ({x) + £(-y) + ((y - u) 

Upon using the definition of $, the right hand side of this equation takes the 
form 



$(x + y; v\ 



<j(y — x)o{y — y)o{x + y) 
a(v-x- y)a(u)a(x)a(y) 



C(y-x)+C(x)+C(-y)+C{y-^)\- (15) 
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After noting the two identities p3 



C{x) + t{y) + t(z)-({x + y + z) 



a(x + y)cr(y + z)a(z + x) 



and 



pO) - p(v) 



a(x)a(y)a(z)a(x + y + z) 
a(y - x)a(y + x) 



(16) 



(17) 



a 2 (y)a 2 (x) 

we find ([15]) simplifies to $(x + y; v)[p(x) — p(y)], where p(x) = —('(x) is the 
Weierstrass p-function. Putting these together yields the addition formula 

<5>(x;v) ${y;v) 

$(x + y;v) = 1 ' > f^—, (18) 

W - fPW 

and consequently a solution of ([!]) with the stated form. Further, from ([18]) we 
see the solution to (|5]) mentioned in the introduction. 

The general solution (fj) involves the two nonzero constants V\ and z/ 2 . Let 
us see how our group of symmetries enables <pi(x) = $(a:; ui) to occur as a limit 
V2 0. Consider the Q action on the general solution <fii(x) of theorem 1 given 
Ax) where 



by (f)i(x) - 
M x ) 



Mx] 
Mx] 



M X ) 
05 (x) 



-v 2 



M X ), 



z/i) 



Now 



and 



lim 0! [x) = $(x; z/^ 



Mx) My) 

M x ) My) 

Thus (^]) arises as the z/ 2 — » of (0). 



lim 

V2~>0 



p(x) - p{y). 



4 Proof of Theorem 1 

Our proof of Theorem 1 proceeds in two stages. First we will use the symmetry 
(]T0| ) to transform (|IJ) into a particularly simple canonical form. This form may 
be immediately integrated to yield a functional equation studied by Bruschi and 
Calogero [[|; by appealing to their result our Theorem 1 will follow. The first 
stage of this process is entirely algorithmic and consequently we may readily 
identify the parameters that appear in our solution. We begin with 
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Definition 1 A point xq G C is said to be generic for ([![) if 

1) 4>k(%) is regular at x for k = 2 . . .5, 

2) <f>i(x) is regular at 2xq, 



3) 



02 Oo) 02(^0) 
0' 3 (^o) 



7^0 



^0. 



04(^o) 04(^0) 
05 (X ) 05 (X ) 

Now let Xo be a generic point. Using at first the matrices U and V of 
transformation ( |10D we may choose linear combinations of <pk (k : 2 . . . 5) such 
that flip becomes 

02 (x) 02 (y) 

03 (X) 03 (y) 



0i(x + y) 



4 (x) 4 (|/) 

0s(^) 0s(y) 



and (for k — 1,2) 



2fc (O) = 2fc+ i(O) = 0, &(0) = 2fe+ i(O) 



(19) 



The arguments of the functions have been shifted to be centred on xo (or 2xo in 
the case of 0i(x)). Here we have set 



and (for k — 1, 2) 

( _02fc(x) 



01 (x) = C0 1 (x + 2x ) 



02fc( X o) 02fc(x O ) \ V 02fe( X + ;r o) 

02fc+l( X o) 02fc+l(x O ) / V02fc+l( X + X Oy 



The constant c here is just the ratio of the appropriate determinants specified 
in ([TOP . We next observe 

Lemma 1 For k = 1,2 we may write 

( j>2k(x) \ = 1 / £ k (x) 



where 7fc(x),^(x) are regular at and 

6(0) = 0, eUO) = 7fc(0) 



1. 



Further, upon writing £&(x) = e AfcX 6(x) = — (f>2 k (0)/2 the function £fc(x) , 

regular at 0, satisfies 



4(o) = g(o) = 0, g(o) = 1. 
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Proof. Upon differentiating £, k (x) = 7fc(x)</>2fc(£) and comparing with £' k (x) 
r Yk{x)<p2k+i(x) we see that 



7ft 0) 4>2k+i(x) -4>' 2k (x) 



lk{x) 



4>2k(x) 



(20) 



The only issue is whether the righthand side of this differential equation is regular 
at x = 0. Using (|T9|) and l'Hopital's rule we find 



y*(o) &+i(o) - &(o) 



7ft (0) 



02ft (0) 



-&(o). 



and so 7ft (x) and hence £ k (x) are regular at 0. Now 7ft(x) is determined by (|20| ) 
given an initial condition which we choose to be 7ft (0) = 1. The remaining initial 
conditions for £ k ( x ) follow from (0). Indeed from 



02ft+l( X ) 



m*hk(x) - ? k (xH(x) 

7ft (x) 



we also find 



(0) = -&(0). 



Now upon writing = e kX £,k(x) with Aft 

statement of the lemma. 



^2ft(0)/2 we obtain the final 

□ 



Thus far we havent used the exponential part of the symmetry fllQ|). Utilising 



this symmetry we set £o(x) = e 



>!{x , 7(xl 



3 2 ( Al - A2 ^7 2 (x)/7i(x). This 



scaling has (upon noting 2X k = 7 fc (0)) the effect of making 7'(0) = 0. Thus we 
obtain 



Corollary 1 At any generic point we may rewrite ([][) using the symmetry 
as 



£a{x + y) = i{x)i{y) 



where for k = 1, 2 

lft(0) = Ift'(O) = 7'(0) 



0. 



fj(g) g(y) 
§(*) §(y) 



lft(0) = 7(0) = 1. 



(21) 



:22i 



Given the complexity of the differential equation 
(|2T| ) simplifies much further. In fact we find 



one may wonder whether 
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Lemma 2 The functional equation and consequently may be written 



as 



d f ji(x + y) \_ d ( i^c + y) 



where d = d x — d y . Further 



6(g)6(y) 
6(g) 



6(g)6(y) 

6(g) 
6(g)' 



(23) 



(24) 



Proof. Taking the logarithmic derivative of (|21f) with respect to 9 = d x — d y we 
obtain 

i{x) i{ y ) 82 (6(^)6(2/)) 92 (6(g)6(y)^ 



7(g) 7(1/) 9(6(x)ei(y) 
Now employing (p2|) one finds 



5 



6(g)6(y)) ' 



(25) 



^ (6(g)6(y) 

and similarly 



|y=0 



&(*)&(v) - f*(*)S(v) 



"6(g) 



d 2 U k ^My))\y=o = -2e k (x). 



Upon setting y = in (P5|) and with these simplifications we obtain the differ- 
ential equation 







7 '(x) + 2£(x) 2&(x) 



with solution 



7(g) 6(g) 6(g) 

7(g) 



c 6 2 (g) 
6 2 (g) 



(26) 



Again using l'Hopital's rule and (^) we find the constant c = 1. Therefore fl2l| ) 
may be rewritten as 



6(g + y) 



Letting y — > we find 







6 




6(y) 




»g(y) 


g 








» 6 2 (y) 


6 




6(y) 










§(v) 



'— ) 



1 V 



1 V 



d{ 



62 (2)6(2/) ' 



(27) 



6(g) 



6(g) 
6(g) 
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as required. Utilizing fl2"3|) we may immediately rewrite ( f2~T|) in the stated form 

O) ^ 

We observe that at this stage the symmetry ([H]) has enabled us to transform 
([J) into the form specified by theorem 1. The solution will follow once we show 
= v k \ Now p3| ) may be immediately integrated to give 

f l(g t y) = U ^ V) + Q(x + y). 

£i(xMy) 6(*)6(v) 

Upon setting = £ 2 {x)/£i(x) and ^(x) = Q(x) / £, 2 (x) this may be rearranged 
into the form 

a l x+ y\ = i + Ux)UvW* + y) (28) 

which is the functional equation studied by Bruschi and Calogero ||. Calling 
upon the general analytic solution obtained by these authors, together with our 
initial conditions (|22|), we findP] l/^ k {.x) = <&(x] v k ) as required. 



5 Proof of Theorem 2 

It is useful at the outset to gather together the various transformations intro- 
duced in the last section: 



4>2k(x) \ f 4>' 2k (x ) <p2k(x ) \ V 4>2k(x + X ) 

hk+l{x) ) V02fc+l( X o) 02fc+lOo) / \<f)2k+l{x + X ) 

1 fUx) \ _ eXkX ( i o\f£ k (x 



7*0*0 \€' k (x)J lk (x)\X k lj\t k {x) 
e XkX ( 1 \ / $(x; v k ) 



lk {x)$ 2 {x-v k )\\k -\)\®{x\v k 



Let us introduce the function 



(29) 
(30) 
(31) 



1 ' 0i(2x o ) <&{x;v 2 ) 1 ; 



/(x) = lk {x)$ 2 {x-v k )e- 2XkX . (33) 

3 For example, from Q we obtain £2 s ) = Ae cx a(ax\uj , to') / a(ax + v\lu , lu') . Using the prop- 
erty <j(ax\au, aw') = acr(ax|w, u>') and the definition of <fr(x; v) we may rewrite this as £2(2^) = 
{A/cr{v/a))e^ c -^ v ' a ^ x I${z\-v/a). Now the x -► limit shows {A/a{u/a))e^-^ v l a ^ x = f. 
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Observe that (p^ ) entails the function f(x) is independent of k, 

7i(x)$ 2 (a;; v^ 2 ^ = l2 (x)&{x- v k )e~ 2X * x . 
With this definition we may rewrite (^9|) and (|3l"l) to give 

/ 2 fc(^ + a;o) 
V 02*+i(z + £o) 



e XkX ( 0' 2 fc(^o) <M^o) \, 

^2k+l 



f(x) V ^fc+l^o) 02fc+l(^ O ) / \ A fc 







/(*) 



z/ fc ) 02fc(^o) 
$'(a;; z/ fc ) (f>' 2k (x ) + \ k (p2k(x Q ) 

$'(x; v k ) (f> 2k+1 (x ) + ^k(p2k+i(x ) 



v k ) 
&(x; v k ) 

(34) 

\ 



and 



u k ) 



f(x)e 



1 



02*0*0) 4>2k(x ) 

02fc+i(^o) 02fc+i(a;o) 



/ 



/(x)e 



V 02*(£ + £o) 
\02fe+i(a3 + x ) 
(35) 



02fc(a; + a;o) 02fc(^o) 
b 2k+ i(x + x ) 4>2k+i(xo) 

X + X ) 2 fc( X o) + Afe0 2 fe(Zo) 

4>2k+l{x + X ) 4>2 k+1 {Xo) + \k4>2k+l(x ) 



\ 



>2k\ 



(j)' 2k (Xo) <l>2k(Xo) 
02fc+lOo) <l)2k+l(Xo) \ <t>2k+l(x + X ) (j>' 2k+1 (x ) + A fc 02fe+l(^o) / 

Now ([32D and (^4|) are of the form stated in theorem 2 provided we can show 
f{x), defined in (p3|), can also be put into the form of the theorem. To see this 
note that (E9|) shows 



hk{x) 



02fc+l(^o)02fcp£ + X ) - 4>2k(Xo)(j)2k+l(x + X ) 
(f)2k+l(x )(j)' 2k (xo) - (p2k(x )(f)' 2k+1 (xo) 



4>2k{x + X ) 4>2k{x ) 
4>2k+l{x + X ) 4>2k+l(x ) 



(36) 



4>2k( X o) <l>2k(Xo) 
02fe+l ( X o) 02fe+l(^o) 



while from ([Jl|) we see 



7* 0) 



$(x; 



2* (a? J 



(37) 



Combining these thus shows 
/(*) 



02fc( x o) 02fc(a;o) 

fak+lM <fi2k+l(x ) 



4>2k{x + X ) 4>2k(x ) 
(j) 2 k+l(x + X ) 4>2k+l(x ) 



(3* 
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as required. Also from fl36|) and the definition X k = — 02 fc (O)/2 we find 



-2A» 



02fe+l(go)j2fc(^o) ~ 02fc(^o)02fc+l(^o) 

2 fc+i(xo)0 2jfc (x o ) - 02fc(xo)0 2fc+ i(x o ) 

(f)' 2k (x + X ) 02fc(x O ) 
>' 2k+1 (x + X ) 02fc+l(x O ) a ._ () 



X In 



(39) 
(40) 



At this stage then we see that if we can determine $(x; v k ) a U of the terms 
in (^9|-[32|) are determined and we obtain the stated expressions for 0i(x), 2 (x), 
03 (x), 04 (x), 05 (a:) and /(x) given in theorem 2. It remains therefore to deter- 
mine the parameters #2, #3 specifying the elliptic functions $(x; i/&) as well as 
ui, vi- To this end we utilise (p5|) to give 



g^Xj Z^fc) 

$(x; i/ fc ) 



02fe(x + X O ) 


02fc(^o) 


02fc+l(x + X ) 


0^+1(^0) 


02fc(x + X O ) 


02fc(x O ) 


02fc+l(x + X ) 


02fc+l(x O ) 



<9<, In 



02fc(x + X O ) 02fc(Z/ + X O ) 

02/t+i (x + x ) 02^+1 (y + x ) 



2/=0 



Upon using (|TJ) to simplify the left-hand side of this equality we thus obtain 
the first equality of theorem 2, 



d y In 



02fc(x + X O ) 02fc(2/ + X O ) 

02/t+i (x + x ) 02fe+i (y + x ) 



C(vk)-((x)-((i/ k -x)-\ k , (41) 



y=o 



and consequently 



d x d y In 



02fc(x + X O ) 02fe(2/ + X O ) 

02fc+l(x + X O ) 02fc+l(2/ + X O ) 



p(x) - p(v k -x). (42) 



y=0 



In fact we have the more general result 



d x d y In 



02fc(x + X o ) 02fc(Z/ + X O ) 

02fc+l(x + X O ) 02fc+l(2/ + X O ) 



d x d y In 



$(x; u k ) $(y; i/ fc ) 



In 



$(x + y; u k ) [ p(x) - p(y) 

p'{x)p'(y) 



p(x + y)- p{v k -x-y) + 



p(x) - p{y)^j 



from which (fE3) arises as the y — > limit. 
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It remains to how that the Laurent series of (fllj) (or equivalently of (|42|)) 
determines the parameters of $(x; v k \ Set 

1 X l+1 

c(u k ) - ax) - au k - x) - \ k = — \ k + Y< Fi jj— w (43) 

X 1=0 V + J- 

or equivalently 

p(x)-p(v k -x) = -^ + J2jX l - (44) 

1=0 

While the coefficients Fi in these expansions depend on k = 1, 2 we will avoid 
including this in our notation: certainly the combinations of these coefficients 
giving g 2 and g% are independent of k. Now the left-hand side of ([43]) has the 
expansion 

Afc - p(v k ) x + p'(v k ) — + (2c 2 - p"{v k ))— + ... 

x 2 3! 



while that of ( |44|) begins 
1 



— + c 2 x 2 + c 3 x 4 + 

ar 



- |pK) - ^p'(^fe) + yp"K) + • • • } 
From either of these we see 

F = -p{u k ), Fx = p'(y k ), F 2 = 2c 2 - p"{v k ), 
whereupon utilising (|67|) we obtain 



c 2 = ^|^ = | and g 3 = 6F*-F?+ 5 -F F 2 . (45) 

Thus, as stated in Theorem 2, we may obtain the parameters of the elliptic 
functions from the Laurent expansion ( f4~3"| ) for either choice of k, the combinations 
of the coefficients in (f45|) being independent of k. The constant terms in the two 
expansions the determine vi, v 2 via F = — p{y k \ g 2 , g^). 

We now have now established all of Theorem 2. It is perhaps useful to 
conclude the section with a lemma that implements the theorem. 



Lemma 3 Let 



d x d y In 



4>2k(x + x ) 4> 2k (y + x ) 
02fc+i (x + xq) <p 2k+1 (y + x ) 



x l+1 



h'(0)h'(x) 
(h{x)-h{0)) 2 



13 



where h(x) 



F 



4>2k(x + X ) 



. Set hi 



h x 1 
h 2 hi 



Fi 



h( k+1 \0) 
(k + l)\h'(0) 



. Then 



hi 1 
h 2 h 1 
h 3 h 2 hi 



hi 


1 








h 2 




1 





h 3 


h 2 


hi 


1 


hi 


h 3 


h 2 


hi 



6 Examples 

We shall now consider the classical addition theorems of the Jacobi elliptic func- 
tions and then the functional equations (|j) and (|7|) as examples of our theory. 
We have collected several standard results pertaining to elliptic functions that 
are of use in our computations in Appendix A. 

Example 1. As a first application of our theory we consider the addition the- 
orems for the Jacobi elliptic functions dn(x), cn(x) and sn(x) where dn(a;) = 
dn(x\m) and so on. These may be cast in the form of ([!]) as 



dn(x + y) 



and 



cn.\x) 


ca'(y) 


cn(rc) 


cn(y) 


sn'(a;) 


sn'(y) 


sn(x) 


sn(y) 


sn 


'x + y) = 



(Jacobi), cn(x + y) 





dn'(x) 


dn'{y) 


1 


dn(x) 


dn(y) 


k 2 


sn'(x) 


sn'(y) 




sn(x) 


sn(y) 



1 1 

sn 2 (x) sn 2 (y) 



sn'(x) 
sn(x) 



sn'(y) 
sn(y) 



(Cayley). 



Let us now apply our theorem to the first equality. The first step is to choose 
an appropriate generic point Xq. This means we wish xq to be a regular point 
for cn(x), sn(x) and dn(2x) as well as 



^ 
^ 



cn'(a;o) cn"(x ) 

cn(xo) cn'(x ) 

sn'(xo) sn"(x ) 

sn(x ) sn'fxo) 



1 — m + m cn 4 (x ), 
1 — m sn 4 (xo). 



Thus we can take Xq = for this example. 
Using (|39D we find 



— 2Ai = d x lncn'fx) 



!x=0 



and 



2Ao = d T lnsn'fx) 



\x=0 



0. 
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Further, with h(x) = fafe) / fofe) = cn'(x)/ cn(x) we obtain 

1 — m + m cn 4 (x) 1 1 — 2m 1 + 14m — 14m 2 



F(x) 



sn 2 (x) dn 2 (x) 



+ 



15 



x 2 + 



while with h(x) = <fi^(x) / <fi§(x) = sn'(x)/ sn(x) we obtain 



F(x) 



m sn ix 



sn 2 (x) 



1 + m 1 



16m + m 2 



x- 
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-x z + ... 



In both cases we find 



92 = -(1 - m + m 2 ) and # 3 = — (m - 2) (2m - l)(m + 1), 
(the required equality providing a nontrivial check) which means 



m 



2m- 1 



ei 



e 2 



and e 3 



-1 



m 



Further, 



2m - 1 



and p(i/ 2 ) 



-1 



m 



3 v 7 3 

Comparison with (|TDD and ( ff2| ) shows co> = K(m), uj' = iK'(m), v\ = K(m) + 
iK'{m) and z/ 2 = K'(m). We may also calculate /(a;) = l/sn 2 (x) and upon 
using ( |76|) our identity may be rewritten as 

$(x + y; K(m) + iK'{m)) 



dn(x + y) 



$(x + ?/; iK'(m)) 

$(x; K(m) + iK'(m)) $(y; iiT(m) + iK'(m)) 

$(x; if (m) + iK'(m))' $(y; if (m) + iK'(m))' 

$(x;iK'(m)) Q(y;iK'(m)) 

$(x;iK'(m)y <S>{y;iK'(m))' 



1 


cn' 


(x) 


cn'(y) 


sn 2 (:r) 


cm 


z) 


cn(y) 


1 


sn' 


(x) 


sn'(y) 


sn 2 (x) 


sni 


x) 


sn(y) 



The second identity may be treated in the same manner yielding v\ = K(m) 
and V2 = K'(m). The third identity is a little different. It may be rewritten as 



$(x + y; iK'(m)) 



1 



sn 2 (x) 



sn'(x) sn'(?/) 
sn(x) sn(y) 



sn(x + y) 



sn 2 (x) 



1 1 

sn 2 (x) sn 2 (?/) 
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$(x;iK'(m)) $(y;iK'(m)) 

$(x\iK'(m))' $(y;iK'(m))' 

& 2 (x;iK'(m)) & 2 (y;iK'(m)) 

1 1 



Now 



$ 2 (x;iK'(m)) - <& 2 {y]iK'(m)) = p(x) - p{y) 

and the required identity follows from the general solution by the limiting pro- 
cedure described in section 3. 

Example 2. We shall now determine the general analytic solution of 



M x + y) = M x )My) + My)M x ) 



M x ) My) 
M x ) My) 



M x ) My) 
M x ) My) 



x), 03 (x) = (f>l(x). The particular case <j>i(x) = 4 (x) was 



where 2 (x) = 
treated in 0. 

Suppose xq is a generic point. Then from 



0^ 



M x o) M x o) 
M x o) M x o) 



2(f)4(x Q )M x o) 



M x o) M x o) 
M x o) M x o) 



we see 4 (x o ) ^ 0, 0s(x o ) ^ and 0i(2x o ) = 20 4 (xo)05(x o ) ^ 0. Further, from 
J), we find 

1 



Ai — A2 



M x o) M x o) 



(46) 



^M x o) M x o). 

Our strategy is as follows. We will first determine u±, u 2 , Xi, A 2 , the pa- 
rameters describing the elliptic functions and the ratio 4 (x + x )(f) 5 (x )/4> 5 (x + 
x o)4>i( x o)- Then from 



>i(x + 2x c 



4 (x + X )M X o) + <pi{xo)(i) b (x + Xq) 

5 (x + Xo)04(Xo) 
4 (x + Xo)4> 5 (x ) 



4 (x + Xo)0 5 (^o) 1 + 



$(x; i/ 2 ) 



we will obtain 

4 (x + Xq) 



20 4 (x o )e 



(A2— Ai)i 



$(x;z/i) 



1 + 05 (X + X O )0 4 (x O )/0 4 (x + Xo)05(Xo) lAj) 

with 05 (x + x ) immediately following. 



(47) 
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Now from (153) we obtain 



1 + — (4h 



where 



02fc(x o ) 02/t+i(a;o) " 02fc+i(x o ) " $(x;z/ fc ) 

Here TVj = 2N 2 and by our assumption that Xq was a generic point these are non 
vanishing. Further, from <p 2 (x) = 4>\(x) and 3 (x) = </>f(x), we see that 

1 + N 1 /D 1 = (1 + N 2 /D 2 f . 

Expanding this shows D\ = (D 2 + N 2 /2)Di which upon using fl46|) yields 

0' 5 (X O ) | = M(^) | | iV 2 (4Q) 

0^o) +A iV 2 $'(x; 
5 (a;o) 2 2 $(x;i/i) 



Suppose that v 2 is finite. Comparing the pole behaviour of each side of ( [SO]) shows 
that i/j = v 2 and consequently that N 2 = 0, a contradiction. The remaining 
possibility is that i/ 2 is infinite which we now show to be a consistent solution. 
This can only happen if the elliptic function degenerates into a hyperbolic or 
trigonometric function and without loss of generality we choose the former. In 
this case 

*(z;u)= ^M^Z^ ^c otw and $ ( ) = _* (50) 
smh kv smh kx smh kx 

Let us the suppose v 2 = oo. Utilising ( [77]) and (|78|) we then must solve 
+ X 2 + k coth KX I = — ^ — - + A2 H V k coth 



^fao) " / V^5(^o) " 2 

< ^ >5 ^ :ro | + A 2 + k coth k(z/x — x) + k coth kx — k coth kz/x 

>5(^o) 2 

This holds provided 

2 4/t 2 <// 5 (x ) iV 2 

iv 2 = -T-; , -r - ; — r + A 2 + — + « coth KV\ = 

smh/tz/x <P5{Xq) 2 
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which determines V\, X 2 and (via Ai in terms of 04(xo), 05(xo), ^4(^0) 

and 4>' 5 (xo). The choice of sign in taking the square-root here is arbitrary (just 
defining v{) and we will take N 2 = —2k/ sinh/«i/i. Substituting these into 
we find 

4>i{x + Xq) <P^{xq) K COth KX — K COth KUi + N 2 /2 

5 (x + Xq) 4>a(Xo) ft coth — K COtllftZ/! — N 2 /2 

= coth(ftr/ 1 /2) tanh«(z/i/2 — x). 



Now employing fl5"0|) shows 



x (x + 2x ) = e (A 2 -A 1+ .coth Wl )^ sinh x ) l(2 x o ) (51) 

smh kv\ 

where the exponential may be rewritten to yield 

\ \ 4-1, ftiM N 2 fao) . ,1 / /o n 

A2 — Ai + KCOthftZ/x = h KCOthKZ/i = + KCOth{KUi/2) 

05 Oo) . N 2 0' 5 (x o ) , . . 

= — H h k coth kv\ = —— — - + Ktanh(/ti/ 1 /2). 

5 (x o ) 2 05(^o) 



We now have the information needed to determine 04 (x) and 05 (x) via ( 47 ) 
which gives 

04 ( X + Xo ) = sinh e ^ ixo y Mxo)+KmtHKVl/2) ) x , } (52) 

sinh(/tz/i/2) 

05 ( X + Xo ) = COsh e fato)/<hto)+«to**W2))x Hxq) (53) 

cosh(/tz/ 1 /2) 

Assembling this provides 

Lemma 4 T7ie general analytic solution to (fjj en by (|57|j, (|^) and 

where xq is a generic point. 

Example 3. We conclude by showing how our results determine the solutions 
of the functional equation (jfy: 

^x(x + y) = V 2 (x + y)<h(x)(f) 3 (y) + V 3 (x + y)0 4 (x)0 5 (y). 

This equation encompasses as particular cases the equations (|j) (with ^ 2 = 
\1>3 = 1, 02 (x) = 04 (x) and 03 (x) = 05 (x)) and (|6|) (with (0 2 (x) = 03 (x) and 
04 (x) = 05 (x)) which have already been discussed. Because of this we will only 
consider the generic case 02 (x) 7^ A03(x), 04 (x) 7^ 70s (x) and ^ 2 {x) 7^ (x) 
(where A, 7, 5 are constants) rather than these limits. Our first step is to relate 
(0) to ©: 
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Lemma 5 The functions $? m (x) (m = 1,2,3) and 4> n (x) (n 
solution of equation (Q) if and only if 



^ 3 (x + y) 
^ 2 (x + y) 



<P2{ X ) 



h(y) 
My) 



an 



,d 



®i(x + y) 
^2(x + y) 



(f>i{x) 



b 2 (x)(j) 5 (x) 

b 3 (x)(j} 4 {x) 



h(y) 
h(y) 



2, 3, 4, 5) give a 



(54) 



My)My) 
My)My) 



(j) 4 (x) 

05 (x) 



k(y) 
h(y) 



(55) 



Proof. Assume first that the functions \l/ m ,0 n give a solution of equation (|7|). 
Then after interchanging x and y in (|7|) and subtracting the result from (^) we 
obtain equation fl54l). Upon substituting the formula for ^(x + y)/^/ 2 (x + y) 
into (|7|) we arrive the formula (|55|). 

In the other direction, let the functions \I/ m ,0 n now satisfy (^4j) and (|55"D . 
Upon writing the right hand side of (|7|) as 



^2(x + y )Mx)My) + ^3(x + y)Mx)My) = 
* 2 (* + y) U 2 (x)My) + %$$M*)Mv) 



(56) 



and using expression (|54| ) for ^(x + y)/ty 2 (x + y) we find the term in brackets 
in (|56| ) rearranges to give precisely the right hand side of (|55|) ; substituting for 
this then yields (|?J) and therefore the required solution. □ 



We may now apply theorem 1 to show that if \l/ r 
solution of (0) then we must have the ratios 



(x 



[m 



fijx) 



$(x ;y u 2 )' V 2 (x) 



c 2 e 



\ 2 x ®(x;ii 3 ) 
$(x;/i 4 ) 



1, 2, 3) give a 



(57) 



Further, because the denominators of (|54|) and (|55|) are the same, theorem 2 
shows that /14 = fi 2 . Theorem 1 also determines the functions <fi n (x) (n = 
2,3,4,5) up to a Q action. In fact, given three functions ^ m [x) {m = 1,2,3) 
whose ratios satisfy fl5"T|) with /z 4 = fi 2 , this is also sufficient to guarantee there 
are functions <fi n (x) {n = 2,3,4,5) for which (|7|) holds true. To see this let us 
substitute these ratios into equation (|7|) to give 



Cl e Xl{x+y) <S>(x + y; /n) = $(x + y; /J 2 )(f) 2 (x)My) + 
c 2 e X2i - x+y) <&( 



[x 



y;H 3 )(f) 4 (x)(f) 5 (y). 



(5f 



19 



We will have established sufficiency once we have shown how to construct the 
functions n (x).This will be achieved utilizing various properties of the functions 
v). 

Lemma 6 The Baker- Akhiezer functions v) satisfy the equations 

$(x + a;u) = - e (««-")+«'')-C(«))*$( a; v ) - ~ a) (59) 

$(— x; a) 

and 

ce 7(^+?/) $(3- + y; ^ + zy 2 ) = $(;£ + y; ^) $(3;; u 2 ) $(y; I/ 2 ) - 

$(x + ?/;z/ 2 )$(x;z/i)<l>(y;z/ 1 ), (60) 

w/iere c = p(i/ 2 ) - p{vi) and 7 = COi) + C(^) - C(^i + v -i)- 

These follow directly from the definition of v) and properties of the Weier- 
strass sigma function; in particular ( pD| ) is a consequence of the 'three term 
relation' of the sigma function ( f2~3| , 20.53, Ex:5\). 
Upon setting x — > x + a in (|60| ) we obtain 

ce 7(a:+2/+Q!) $(x + y + a; ^1 + v 2 ) = $(x + y + a; z/i)$(x + a; v 2 )®(y, v 2 ) 

-$(x + y + a; z/ 2 )$(x + a; i^)$(y; ^1) 

(61) 

Now by substituting ( |59"D in ([H]) and setting /ii = z/i + v 2 — a, /x 2 = v\ — a 
and /i3 = v 2 — a one obtains after some rearrangement 

(62) 

for appropriate constants d , c", A', A". This is precisely of the desired form fl5"8|) . 
Therefore we have shown 

Theorem 3 Given functions ^ m (x) (m = 1,2,3), there are functions (j) n {x) 
(n = 2,3,4, 5) for which the functional equation (Qj is true if and only if the 
following ratios take place: 

jMg) = r r X lX ®&»l) gg(g) = r W ^O^s) , fi „v 

* 2 (x) 1 $(x;/i 2 )' * 2 (x) 2 $(x;/i 2 )' ^ 
where c m , X m , (m = 1,2) and // n; (n = 1, 2, 3j are free parameters. 
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A Elliptic Functions 

A.l The Weierstrass elliptic functions: 

The Weierstrass elliptic functions are based on a lattice with periods 2u and 2u/, 
where $t(u'/u) > 0. They satisfy the homogeneity relations 

cr(tx\tu>,tu>') = ta(x\uj,uj'), C,{tx\tuj,tuj') =t~ 1 ( > {x\uj,uj'), (64) 
p(tx\tuj,tuj') = t~ 2 p(x\uj, uj'). (65) 

Here ((x) = (\aa(x))' and p(x) = —C,'{x). These homogeneity relations mean 
our function <3>(x; v) = v\u,u>') satisfies 

tv\tu, tu') = t _1 $(x; v\u, uj'). 

The parameters g 2 and g% alternately used to describe the elliptic function 
are given by 

g 2 = 60 IT 4 , ^3 = 140 J2' ^ 

where f2 = 2mw + 2nuo' . 

The Weierstrass p-function, = p(x\uj,uj') = p(x; g 2 , gz), satisfies the 

differential equation 

p' 2 (x) = 4p 3 (x) - # 2 p(:r) -g 3 = A[p\x) - 5c 2 p(x) - 7c 3 ) (66) 

where c 2 = || and c 3 = ||. This means 

p"(x) = 6p 2 (x) - 10c 2 . (67) 
The terms c& (A; > 4) in the Laurent expansion of the Weierstrass p function, 

Ob 

1=2 

are expressible in terms of c 2 and C3. 

If ej = ej(o;, uj') (i = 1, 2, 3) denote the roots of the cubic 

4x 3 - g 2 x -g 3 = 0, (69) 

whence 

p l2 {x) = 4(p(x) - e x ){p{x) - e 2 )(p(x) - e 3 ), 
the half-periods a;, are defined by 

p{uJi) = where uj\ — uj, uj 2 = uj + uj', and uj 3 = uj'. 
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Clearly 

ei(tu, tui') = t- 2 ei(cu, J). (70) 

Assuming that gi and are real, and that the discriminant of (|69"D is positive 
(the case of interest in this paper) the e^'s are real and may be ordered e\ > 
e 2 > e 3 . 



A. 2 The Jacobi elliptic functions: 

The Jacobi elliptic functions are characterised by a parameter m. Thus for 
example sn(x) = sn(x|m) has periods AK(m) and 2iK'(m), where K{m) is the 
complete elliptic function of the first kind. The e^'s are related to the parameter 
m of the Jacobi elliptic functions by 

2 — mK 2 (m) 2m — lK 2 (m) , — l — mK 2 (m) , , 

— y> « =— — y a,id ca =— — <n) 

Thus 



_ 4,. o.KHm) 4 , . ^AT 6 (m 

and 



<? 2 = ^(1 - m + m 2 )^p, g 3 = -Um - 2) (2m - l)(m + l)^p. (72) 

o UJ Z l UJ 



uj _ iK'(m) K{m) 

w K[m) V e i _ e 3 

We find 

i\ a *.t \ cn(ax) , .,. dn(ax) 

$(x;u/) = — - — r, §(x;uj) = a — - — (■ and $(x;u + w ) = a ■ 



sn^axj sn{ax) sn^axj 

Here a = y/ei — converts the periods of $ based on uj, uj' to those of the 
Jacobi functions based on K(m), iK'(m). Equally we may write this as 

$(z; tu/ fcc;, to;') = — —. — r, with t = yfe x - e 3 = v 7 . (74) 
sn{x\m) uj 

Thus, with these periods in mind, we write 

<&(x;iK'(m)) = — t~T) &x]K(m)+iK'(m)) = ^Q (75) 
sn(x) sn(x) 

$(x]K(m)) = ^Ff (76) 
sn a; 
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A. 3 Degenerations: 

When the discriminant of (|69|) vanishes one (or both) of the periods of the elliptic 
function vanishes yielding hyperbolic, trigonometric (or rational) functions. If 

e i = e 2 = c, e 3 = —2c (and so g 2 = 12c 2 , g 3 = —8c 3 ) we then have 

a(x; 12c 2 , -8c 3 ) = e" K x /6 and p(x; 12c 2 , -8c 3 ) = — + 



K 3 sinh 2 kx 

where k = \/3c. In this case 

$ = K smh K (v-x) eXKCOthKV = k ^ _ cQth &XKCOthKU 

smh ki/ sinh k,x 

(77) 

$'(x; z/) = — k, $(x; i/) (coth k{v — x) + coth kx — coth /tz/) . 



In particular 



k $ (x; oo) 

®(x; oo) = — , and — = — k cothra. (76 

sinh ki $(x; oo) 
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